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Abstract — In this article, we study the discrete time 
Geo/G/l queue with negative customers and multiple 
working vacations. Under the RCH disciple, the arrival of a 
negative customer does not receive service and removes the 
positive customer who is being served. Using embedded 
Markov chain and matrix analytic approach, we derive the 
probability generating function (PGF) and mean of the 
stationary queue length at the departure epoch. From the 
process of the proof and the results, we also obtain the 
probabilities that the server is in idle, working vacation, and 
regular busy period, respectively. Finally, we discuss 
numerical results. 

Keywords — Multiple vacation, Negative customers, Discrete 

time, PGF 

I. Introduction 

In the past twenty years, discrete-time vacation queue model has 
received a rapid increase of research attention both in queuing 
theorists and practitioners because of its wide applicability in the 
study of production and inventory systems, the performance 
analysis of digital communication systems, and 
telecommunication networks including the Broadband Integrated 
Services Digital Network (B-ISDN), Asynchronous transfer 
mode (ATM), and related technologies. 

Queues with negative customers have been studied over 
decades, and the work about negative customers without working 
vacation in discrete time can be found in Atencia and Moreno[l, 
2], where the authors considered the single server discrete time 
queue with negative arrivals and various killing disciplines 
caused by the negative customers. And Jinting Wang and Peng 
Zhang[3] discussed a discrete time retrial queue with negative 
customers and unreliable server. Park and Yang [4] analyzed the 
Geo/G/l queue with negative customers and disasters. 
However the literatures are still much fewer in the model of the 
Geo/G/l queue with negative customers. 

For discrete-time queues with vacation policies, Zhang and 
Tian [5] treated the discrete time Geo/G/l system with variant 
vacations in 2001. In this system, after serving all customers the 
server will take a random maximum number of vacations before 
returning to the service. Meanwhile, the discrete time queues 
exactly analysis with different vacation policies can be found in 
Takagi [6], Tian and Alfa [7], Zhang [8], and references therein. 
In these papers, they all assumed that the server stop service 
completely during the vacations. Recently, Li [9] studied the 
discrete time Geo/G/l working vacation queue and its 
application to network scheduling. But there is no negative 
customer to join discussion. The negative customers are the equal 
of all interference elements in real life. For example, in the bank 
queuing service process if some customers have to stop queuing 
due to special reasons of the individual, then we can seen these 



personal special reasons as negative customers; Another 
example, in the process of cell transfer if the interference of a 
virus causes information element transmission failure, then some 
kinds of virus can be seen as negative customers. Hence vacation 
queue with negative customers is relatively reasonable in real 
life. 

On base of analysis above, in this paper we deal with the 
discrete-time Geo/G/l queue with negative customers and 
multiple working vacations. Under the RCH disciple, the arrival 
of a negative customer does not receive service and removes the 
positive customer who is being served. Yet if the system is empty 
and the arrival of the negative customers, the negative customers 
will vanish and have no impact on the server. Once a positive 
customer arrives and the server is idle, the service of the arriving 
customer commences immediately. If the server completes all 
customers' service and finds the system becoming empty, it will 
leave for a vacation. During the vacation period, if there are new 
positive customers arriving, the server will provide service for 
them by their arrival order at a slower rate. At a vacation 
completion instant, if the server finds customers in the system, he 
will come back to the normal working level and take the service 
with the normal service time. Otherwise, it immediately proceeds 
for another vacation and continues in this manner until he finds 
positive customers in the system. If one customer is still being 
served when the vacation ends, we will assume that this customer 
will begin to receive the service as a new one. 

The organization of the paper is as follows. In Section 2 we 
give a model description and embedded Markov chain. In Section 
3 we use embedded Markov chain and matrix analytic approach 
and derive the probability generating function (PGF) and mean of 
the stationary queue length at the departure epoch. From the 
process of the proof and the results, we also obtain the 
probabilities that the server is in idle, working vacation, and 
regular busy period, respectively. Finally, we give a special case. 
In the Section 4, we discuss numerical results. 

II. Model description and embedded Markov chain 
The features of the discrete time Geo/G/l queue with negative 
customers and multiple working vacations studied here are as 
follows: 

(1) We consider a discrete-time queuing system where the time 
axis is segmented into a sequence of equal time intervals. It is 
assumed that all queuing activities (arrivals and departures) occur 
at the slot boundaries, and therefore they may occur at the same 
time. For mathematical clarity, we suppose that a potential 

customer's arrival occurs in(«, « + ), n =0,1, 2,L , and a potential 
customer's departure takes place in , n = l,2,L , where 

n = lim(n-A«) and« + = lim(« + A«) . This model is known 
as the arrival first system. 



IJER@2013 



Page 393 



International Journal of Engineering Research (ISSN : 2319-6890) 

Volume No.2, Issue No.6, pp : 393- 398 1 Oct 2013 

(2) Two types of customers, positive and negative, arrive working vacation model, any service completion may occur 
according to geometrical arrival processes with the parameters during a service period or a working vacation period. Thus, we 
p and q , respectively. Here we have to state the assumption define a process 

regarding the order of these concurrent events because a positive X n = (L n , J n ) , n > 1 , 

customer's arrival and a negative customer's arrival can where j symbo i izes the system state at the nth customer 

simultaneously occur at the same slot boundary. We also assume , , . „ , 

. departure epoch ( / =0 it the system stays in a working 
that the potential negative customer s arrival occurs at t = n 

before the potential positive customer's arrival. Under the RCH vacation P eriod ' or J << =1 if the s y stem sta y s in a bus y service 

disciple, the arrival of a negative customer does not receive period). Then [X n ,n>l} is a two-dimensional embedded 

service and removes the positive customer who is being served. Markov chain of our queue sys tem and its state space is 

However, if the system is empty and the negative customers _ ,„ „,„,,, N , , . „ 
arrive, the negative customers will vanish and have no impact on v ' t v J ' J > 

the server. The server will stay in the vacation period if there are no 

(3) Once a positive customer arrives and the server is idle, the customers in the system. 

service of the arriving customer commences immediately. Let a } , j > be the probability that there are j customer 

Otherwise, the arriving customer joins the waiting line to be . , , . c T . c e 

b J b arrivals during the service time S b . Let the service time S b of a 

served. Let customers be served at slot boundaries \n,n ), .. . , , . , . . , . , 

V ) positive customer be equal to k in the regular busy period. Then 

n = 0,1,2,L .Further, the service of a customer takes an integer we get the following two cases: 
number of slots, which implies that customers also leave the (i) The service of a positive customer ends; 
system at slot boundaries, and so suppose departure the sever (ii) The service of a positive customer doesn't end, but he be 
m{(n + \)-,(n + l)\, removed by a negative customer. 

\ v ' v '/ Then we can get that 



n = 0,1,2,L .Every customer's service time S b is independent and _ ^ wp| H\- \ k ~' 



e sn . y(i-/>r (i-<o + e ft™ . p j (i-?r(i-?) 9 

obeys the general discrete distribution during regular busy period. 4=™a,j) \JJ i=™xa.j) \JJ 

Its general distribution, probability generating function (PGF), \ m(^\'( \ k ^ J ( \ k 

andmeanare^,fc=l,2,L ,G b (z), and E(S b ), respectively. =( 1 + ^^2 ?i [jJP i l ~P> K 1 ^) • J^. 

(4) If the server completes all customers' services and finds the And the sum of all possible arrival probability in the service 
system being empty, it will leave for a vacation. During the time S b is 

vacation period if there are new customers arriving, the server „ „ « /» n 

provides service for them by their arrival order at a slower rate. ^a =^(l + q) g ( k ' } p' (l-/?)* j (l-^)* 

Service time S v of every customer is independent and obeys the j=o j=o *=max(ij) \J J 

general discrete distribution in working vacations. Its general =(l + q)G b (q). 

distribution, probability generating function (PGF), and mean 

(2i , ,~ T „ i \ . „ / a \ ., Since the sum of all possible arrival probability during 

are g h , k =1,2,L , G l ,{z), and E[S V ), respectively. r 1 

(5) At a vacation completion instant, if the server finds the service time S b is equal to 1, one has G b(l) = ^r~ • 
customers in the system, it will come back to the normal working ^ 
level and take the service with the normal service time. On the Hence their PGF and means are 

other hand, if it finds no customer, it immediately proceeds for A( z ) = f> ; V =(l+4)EE^f*]W i l ~ p)"'' =(l+q)G b ( Pq+ pqz) 
another vacation and continues in this manner until it finds x> j-o*-j w 

customers in the system. If one customer is still being served , 

when the vacation ends, we will assume that this customer will (l + q)qp ,- -j. 

begin to receive the service as a new one. ^'( z )|„i = = P> ^"( z )| z= i = i} + l)\1P) G h yqy 

(6) There is only a server in the system, and service rule is first ^ b 

come first service, and the server can serve only one customer at Similarly, let b.{j = 0,1,L ) be the probability that V > S v and j 

a time. The vacation time V obeys geometrical distribution with customers arrive during the service time S„ and v, (j = 0,1,L ) be 

the parameter . . . , 

^ ^ the probability that V < S v and j customers arrive in the 

p(y=k) = 99 k \k = \,2,h ,v(z) = -—^->E(v) = - -j- vacation time V . Therefore, we can denote 

i-fz (1-6') , x 

We assume that inter arrival times, the normal service time, bj =(l + ^) ^ gf } \ \p J (l-p) ~' (\ — q) (l-^) , 



and working vacation time and normal vacation time are t=max(ij> \J ' , 

mutually independent. Thereinafter, for any real number » » 

, e [0,l],wedenotex=l-x. ^ =( 1 + ^) S Z^(l^) . P 1 {l~ pf i^V ■ 

Let L n be the number of the customers at the instant of the And so we have that 

nth service completion or customer departure, and J the state of , . / — \ , ,/ /-\ / — \\ 

Z^ = ( 1 + 9) G »M'Z v ;=( 1 + 9)( G vW- G vM))- 

server both at a slot division point t + , t = 0,1,2,L . For the ;=o j=o V V / // 
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Since the sum of all possible arrival probabilities in the service 

time S v equals to 1, we get that G v (q) = — . Thus, we have the 

K ' l + q 

PGF and means of them in the following: 



= (\ + q)G,\eq(p + pz)), 

B»(z)\^=(l + q)(eqp) G v " (Oq) 

co 



CO GO 



sn(k 



=i(i+9) x x^wr . (/«y(i-/»r(i-frf 

7=0 jfe=max(l, _/) n=A \J J 

l-0lp + pZ) L J 

\-9[p + pz) 

j 

Let c. =^]v l fl._ i , 7 ^0. Then c ; . represents the probability 

/fc=0 

that V < 5 V and j customers arrive in the vacation time V and 
service time S b , and we have 



oo co J 



co oo 



Xs =ZZw* X>- X" =i-(i+9)g v (^), 

7 =0 7=0fc=0 k=0 j=k 

co co j go oo 

c ( z ) = H c j zJ =XX v * fl y-* zy =X v * z *X fl ;-* z ^ 

j=0 ;=0*=0 t=o j=* 

= A( Z )V( Z ), 

c(C 1 = A '( 1 ) y ( 1 ) +A ( 1 ) y '( 1 ) 

(l + q)G v (0q)]-fi, 



v 0j 



C"(z)| 2 i = A"(l)V(l) + 2A'(l)V'(l) + A(l)V"(l) . 

Now, we consider the one-step transition probabilities 
of (L n , J n ) , « > 1 , and get the following results. 

lfX n =(m, l), m>2,then the remaining customers are those 
m-1 customers and new arriving ones during the service time S b 



at the next departure epoch, and thus X n+1 ={m — \+ j, l) with 
probability a Jt m>2, j > . 

If X n =(l,l), then two cases happen. On the one hand, the 
remaining customers are those that have arrived during the 
normal service time S b at the next departure epoch, on the other 
hand, no customers arrive in the normal service time S b and the 
system becomes empty. Thus, X n+1 =(j',l) with 
probability a t , j > 1 and X n+l = (0, 0) with probability a . 

If X n =(m,0) , m>2, then there also are two cases at the 
next departure epoch. 

Casel: if V > S v , after this departure, the system stay in the 
vacation period and the remaining customers are those m— 1 
customers and new arriving customers in the service time S v . 
Hence 

we have X n+l =(m-l+ j,0) with probability b t ., j>0. 

Case 2: if V < S v , next departure happens in the normal 
service period and the remaining customers are those m— 1 
customers and new arriving customers during the vacation time 
V and the service time S b . By the definition of bj and Cj , we 

have X n+l = (m — 1 + j,l) with probability v . , j > . 

lfX n = (m,0) , m = 0, 1, then by the similar analysis, we get 
that Z n+1 = ( y',0) with probability 

^ , j>\ , X n+1 =(v,l) with probability c j , y>l , 
andX n+1 =(0,0) with probability b + c . 

Using the lexicographical sequence for the states, the transition 
probability matrix of (L n ,/ n ) can be written as the Block-Jacobi 
matrix 

B B { B 2 B 3 ■■■ 
C Aj Aj A 3 ■■■ 

p = Ai A A ■" 
A> A 

where B =b +c ; B, =(fe j ,c f ) , z>l ; C =(b +c ,a ) 1 ; 

c,. 



A = 



,/>o. 



And the matrix P satisfies the following conditions 

B a + £b jC = 1 ; C A,e = e ; £ 4* = e, e = (l, l) . 

Next, we give out a lemma which is crucial to prove the 
theorem in the sequel. 

(l + q)qp 

Lemma 1. If p = - — <1 and 0>O , then the matrix 

co 

equation G = ^ A ( G' has the minimal nonnegative solution 

G= ' ' 
1 
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where y is the unique root in the range < z < 1 of the equation 
z = (l + q)G v {Oq(p + pz)). 



Z = {l + q)G v (eq(p + pz)) 
has a unique root y, < y < 1 , and so r n = y . 



Proof. Since all A i 's are upper triangular, we can assume that Finally, taking r u = y and r 22 = 1 into the second equation, we 
the minimal nonnegative solution G has the same structure as easily have that r 12 = l—y. 

In view of the structure of matrix G , G is a stochastic matrix. 



If i > 1 , we have 
G = 



r,. 



r' 



Substituting G, G' , and A i into the matrix equation G = ^ j A i G' . 

i=0 

we obtain 

1i =Y J b , r L={ ] - + ( l)G v [Oq(p + pr n )^ 



oo i— 1 



i=0 ;'=(> /=() 

r 22 =Xa,r4=(l + 9 )G i ( ? (p + pr 22 )). 

Firstly, we consider the equation 

Z = {l + q)G b (q(p + pz)), 

andte f(z) = (l + q)G b (q(p + pz)). 
Then 

, N / X / — \ G f(^) 

/(0) = (l + f)G fc M ) = -W. 

/(l) = (l + ^)G fe P) = l, 

and so 0</(0)</(l). 

Moreover, for any z , < z, < 1 , we get 

f'(z) = (l + q)pqGi(q(p + pz))>0 
f"(z) = (l + q)(pq) 2 G b " (q(p + pz)) > . 

Since /'(l) = (l + q)pqG h ' {q} = p< 1 , it follows that the 

equation z = /(z)has the minimal nonnegative root z = 1 , and 
thus r 22 = 1 . 

Secondly, letg(z) = z-{l + q)G v [Oqi^p + pz^ . Then we can 
get that 

g(0) = -(l + q)G v (dqp)<0 

G v (e~ q ) 

g{\) = \-{l + q)G v {6q) = \ V=V >0 - 

Further, since 

y = (l + q)G v (0q(p + p Z )) 
is an increasing function for < z < 1 . Thus, the equation 



III. Queue length and main results 
Let (L, J) be the stationary limit of X n = (L n , 7 n ) . Then the 
stationary distribution can be defined by 

=P(L = k,J = j) = Urn P{L =k,J n = j) , (k, j) e Q ; 

^ = (^ ,^ 1 ,^ 2 ,L ), =^oo> ^ ={x k o,x kl ) ■ 

Since the evolution of the chain is governed by the one -step 
transition probabilities given by matrix P , the Kolmogorov 
equations nP = n together with the normalization condition for 
the stationary distributions are given by 

= ( b o + c o ) + K f a = x (b +c ) + k w (b +c ) + 7v n a ,{Y) 



x k =x Q B t + Y j 7r J A k+l _ J ,k>l, 



(2) 



7=1 



(3) 



From Eq. (2), we can obtain the row vector generating function 

go co co k+l 



^o£^z i +-0(z)A*(z)-^A,- 



k=\ z 

co 

where A* (z) = '^ j A l z 

k=0 

Furthermore, 



B(z) C(zj 
A(z) 



z(tf o fl(z)-fc (>c + x m ), >c (C(z)-l) +b {k„ + n K )) 

[»,»(;)-''.(''."■)] | C(z)[^B(z)-t,(^-nt l0 )] | x a (C{z)-\) + b„(n a +x t 



z-B(z) -C(z) 
z-A(z) 

1 



C(z) 



B(z)) (z-A(z))(z-B( z )) 
1 



(z-B(z)) 



(z-A(z))(z-B(z)) 



(z-A(z)) 
(5) 



Now, we give out the main theorem and its proof. 



Theorem 



(l + q)qp 
1. If p = ± l2^L < i t 



p = '—— < 1 , then the stationary queue 

M„ 

length L at the departure epoch has the following generating 
function: 

A(z)(z-l){z-B(z)) + z(y-z){A(z)-B(z)-C(z)) 
L \ z ) = ~ — , t_\\(_ — TTi~Zv\ ^o> 



(z-A(z))(z-B(z)) 



(6) 



where 
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(l-p)(l-(l + q)G v (qe)) 



^~r)f - U+(i-r) 



p+ 



P 9 



J J 



(\ + q)G v (q9) 



(3) The probability that the server is in vacation is 
P(j=0) = x +®(z)\^e 1 

(i-p)(i-r) 



A(z) = {i + q)G h (pq + pqz),B(z) = {\ + q)G v (eq(p + pz)), 

C(z) = A(z)V(z),V(z) = — -£ -Jl-B(z)] . (7) 

l-9[p + pz) 

Proof. The PGF L(z) of the stationary queue length at the 
departure epoch can be derived 
L(z) = 7r +t>(z)e 

z[x Q B(z)-b (x +x l0 )] 

=n * + i^Hz)) 



1+ (W) 



p+ 



p0 
e 



(\ + q)G v (q0) 



v \ ~ J) 

(4) The probability that the server is in regular busy period is 
P(y = l) = ( z )| , 



j v 



P 9 



1- 



P 



\ + {\-r)\p + I f\ (i + ,K; 



Where e, = (l, 0) r , e 2 = (0, l) 7 



zC(z)[w B(z)-6 (>r +*i )] z[^ (c(z)-l) + fc (^ +^ 10 )] Corollary 2. Under hypothesis of Corollary 1, we have the 

(z-A(z))(z-B(z)) + (z^X(z)) mean ° f thC Stati ° nary i qU6Ue l6ngth 1 : 

E(L) = 



7i a \-zA{z) + A{z)B{z) + z 2 B{z)-zA{z)B{z) + z 2 C(z)\ 
: (z-A(z))(z-B(z)) 
+zb (x +x 10 )(A(z)-B(z)-C(z)) 



^-r)f- U+(i-r) 



(1 + 9 )G V ( 9 6») 



(z-A(z))(z-B(z)) ■ (8) 

By Lemma 1, we know that /(0<7<1) is the unique root 
the equation 

Z = (l + q)G v (0q(p + pz)) = B( Z ). 
Hence the denominator of Eq. (8) equals if z = y , and so the 
numerator of Eq. (8) also equals if z = y ■ Substituting z = y 
into the numerator of 
Eq. (8), we can get 

yjv = b ( tt +n w )^>n w = i^hl ^ (9 ) 



{\ + p{y-2) + p 



-\+{y-\)+{y-\)\p+ 



op 



+ (\ + q) 2 (qp)\y-\)G\qe)G:(q) + [\-(\ + q)G v (qe 



p0 



9 



p0 



p-N/' + t-^ + tH-H 







2 



> . 



In this place, we will list a special case of our results which is 
consistent with Li [9]. 

Case: If we consider the case<7 = 0, the model converts into 

the discrete time Geo I Gil queue with multiple working 
vacations, and the stationary queue length L : 
Taking Eq. (9) into Eq. (8), the Eq. (8) can be converted into Eq. 

(6). Using Eq. (8) and the normalizing condition L(l)=l,Eq. (7) A (z)(z-l)(z- B(z))+ z{y- z)(A(z)- B(z)- C(z)) 

follows. L W= (z-A(z)){z-B(z)) ^ 

Corollary 1. Assume that the discrete time Geo I Gil queue 
with negative customers and multiple working vacations is in the wnere 
(l + q)qp 

steady state. If p = < 1 , then we have the following 

Mb /T| , =- 

results: 

(1) The probability that the server is idle is 
(l-p)(l-(l + q)G v (q0)) 



(l-/>)(l-G, 




) 




i+(i-r)( 


p_ 


< 




G, 





P + 



p0 




(l + q)G v (q0) 

J J 



V(z). 







l-0(p + pz) 



(2) The probability that the server is in working vacation is 



{l-p)[{l + q)G v {q~0)- Y \ 



A{z)=G b {p + pz) 
C(z) = A(z)V(z). 



B{z) = G v (e(p + P z)) 



P 



l + (l-y)^- l + (l-y) P + ^ (l + q )G v 



P 



IV. Numerical results 
In this section, we present some numerical results to study the 
effect of the varying parameters on the main performance 
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characteristics of our system. The values of the parameters are 
chosen so as to satisfy 
(\ + q)qp 

p = - — < 1 . And in the analysis above, we obtain the 

expected queue length in the steady state. Hence we will present 
numerical results to explain that our model is reasonable by the 
expected queue length. 



Fig.l . 



— ii=oni 

---- q=0 1 

q=0 15 



D D.1 0.2 0.3 D.4 0.5 0.6 0.7 0.8 0.9 1 
^=0,3,8=0.22 ,y=0.05595 ,^=0.6) 




According to the expression of the expected queue 
length E(L) , we show the effect of the service rate of working 

vacation ju r on the queue length. Figure 1 is two dimension 
contrast figures of the expected queue length 
when q = 0.01, ,q = 0.1, andg =0.15 . Since 



(l + q)qp 

- — < 1 , we get 



get that the effective image is the above 



zero. By the image in Figure 1, the number of the customers in 
steady state decreases along with the increase of /u r . And we also 
find that the number of the customers in steady state also 



decreases along with the increase of the negative customer's 
arrival rate q when the service rate of working vacation ju r is 

constant. In other words, if the parameters p, q, and ju b are fixed 

in this model, then negative customers arrival rate q is more 

higher, the number of the customers decreases more faster under 
the steady state, and while the system is more likely to empty. By 
Figure 2, we find that the number of the customers in steady is 
increasing when the arrival rate of the positive customer p is 

increasing process, which is in accordance with the queuing 
system in real life. 
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